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Abstract. We review the traditional derivation of the fluid-dynamical equations from kinetic theory accord- 
ing to Israel and Stewart. We show that their procedure to close the fluid-dynamical equations of motion is 
not unique. Their approach contains two approximations, the first being the so-called 14-moment approxi- 
mation to truncate the single-particle distribution function. The second consists in the choice of equations 
of motion for the dissipative currents. Israel and Stewart used the second moment of the Boltzmann equa- 
tion, but this is not the only possible choice. In fact, there are infinitely many moments of the Boltzmann 
equation which can serve as equations of motion for the dissipative currents. All resulting equations of 
motion have the same form, but the transport coefficients are different in each case. 
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1 Introduction 

Fluid dynamics is an effective theory to describe the long- 
wavelength, low-frequency dynamics of macroscopic sys- 
tems. In non-relativistic systems, the Navier-Stokes equa- 
tions are able to describe a wide variety of fluids, from 
weakly interacting gases, such as air, to liquids, such as 
water. On the other hand, the theory of relativistic dis- 
sipative fluid dynamics has not yet been completely es- 
tablished and remains a topic of intense investigation. For 
dilute systems, the derivation of fluid dynamics can be in- 
vestigated starting from the relativistic Boltzmann equa- 
tion. 

Chapman-Enskog theory [I] is a well-known approach 
to derive fluid-dynamical equations from the Boltzmann 
equation. In this approach, the single-particle distribution 
function which is the solution of the Boltzmann equation 
is expressed in terms of an expansion in gradients of the 
primary fluid-dynamical variables, i.e., chemical potential, 
temperature, and velocity, each term containing a differ- 
ent power or order of derivatives. This leads to a series 
in powers of the Knudsen number, Kn = X/L, the ratio 
of the mean-free path of the particles, A, and a charac- 
teristic macroscopic length, L. As is well-known, to ze- 
roth order this method leads to the equations of ideal 
fluid dynamics. To first order in Knudsen number one ob- 
tains the Navier-Stokes equations of fluid dynamics. To 
second and higher order in Knudsen number one obtains 
the Burnett and super-Burnett equations. However, rela- 



tivistic Navier-Stokes theory, as well as any higher-order 
truncation of the relativistic Chapman-Enskog expansion 
is unstable and, consequently, unsuitable to describe any 
relativistic fluid existing in nature [2JI31II]. 

The source of such an instability is well understood 
in the relativistic case: it comes from the acausality of 
Navier-Stokes theory 011]. Therefore, a consistent and 
stable theory of relativistic fluid dynamics must also be 
causal. Causal fluid-dynamical equations were first derived 
from kinetic theory by H. Grad [5], for non-relativistic 
systems via the method of moments. In Grad's original 
work, the single-particle distribution function is expanded 
around local equilibrium in terms of a complete set of 
Hermite polynomials [§]. Fluid dynamics is obtained by 
explicitly truncating this expansion, expressing the distri- 
bution function in terms of only 13 moments: the velocity 
field, the temperature, the chemical potential, the heat 
current, and the shear-stress tensor. Due to this trunca- 
tion scheme the method became known as the 13-moment 
approximation. In non-relativistic systems, the correction 
to the equilibrium pressure, the bulk viscous pressure, van- 
ishes and is not included in the 13 variables. 

The generalization of Grad's moment method to rel- 
ativistic systems is non-trivial. One reason is the lack of 
a suitable set of orthogonal polynomials to replace the 
Hermite polynomials 7,8,9 . Despite this problem, rela- 
tivistic generalizations have been given by several authors 
[TUl[TTl[lljlll3lH^] . One of the most well-known works on 
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this topic was done by Israel and Stewart (IS) [8"ll9"UT5] . 
In this approach, the single-particle distribution function 
is expanded in momentum space around its local equi- 
librium value in terms of a series of (reducible) Lorentz- 
tensors formed of particle four-momentum fc M , i.e., l,fc M , 
k^k", ■ ■ ■ . The procedure adopted by Israel and Stewart 
is very similar to Grad's: this expansion is truncated at 
second order in momentum, leaving 14 moments and 14 
coefficients in the distribution function to be identified, 
the so-called 14-moment approximation (in the relativis- 
tic case, the bulk viscous pressure does not vanish, leading 
to one additional moment when compared to Grad's orig- 
inal approach). However, since the expansion is not real- 
ized in terms of an orthogonal set, the coefficients of the 
truncated expansion cannot be immediately determined. 
For this reason, Israel and Stewart chose a set of con- 
straints to express the expansion coefficients in terms of 
the main fluid-dynamical variables. Furthermore, since the 
zeroth and first moments of the Boltzmann equation are 
the usual conservation laws, it seemed natural to choose 
the next (the second) moment of the Boltzmann equation 
to augment and close the conservation equations. 

Nevertheless, this choice is ambiguous, since, once the 
14-moment approximation is applied, any moment of the 
Boltzmann equation will lead to a closed set of equations 
|16U17j . Therefore, inconsistencies may arise because of an 
ambiguity in the choice of the moment equation for clo- 
sure. Recently, it was confirmed that, at least for some 
cases, the IS equations are not in good agreement with 
the numerical solution of the Boltzmann equation 18,19, 
20,21,22,23 . Also, the transport coefficients obtained by 
Israel and Stewart do not coincide with quantum-field the- 
oretical calculations [24] . 

In this paper, we review the derivation of the fluid- 
dynamical equations from the Boltzmann equation using 
the 14-moment approximation, but from a different per- 
spective. First, the single-particle distribution function is 
expanded around equilibrium in terms of an orthogonal 
basis. This allows us to determine the coefficients of the 
expansion without assuming additional constraints. Then, 
we show how the 14-moment approximation emerges from 
such a complete moment expansion. Second, we obtain 
the fluid-dynamical equations for an arbitrary moment of 
the Boltzmann equation and discuss the ambiguity of the 
14-moment approximation. We explicitly show that the 
form of the equations is always the same, regardless of 
the choice of moment of the Boltzmann equation, but the 
values of the transport coefficients are different. Thus, we 
explicitly demonstrate that the traditional 14-moment ap- 
proximation applied to the Boltzmann equation is not able 
to provide a unique set of fluid-dynamical equations of mo- 
tion. In Ref. [T7] it is demonstrated how one can resolve 
this ambiguity by including higher moments of the single- 
particle distribution function, but this is not the subject 
of this paper. 

The ambiguity in the 14-moment approximation is ex- 
plicitly demonstrated by calculating the transport coeffi- 
cients in the 14-moment approximation for a classical gas 
of massless particles with a constant cross section. We use 



two different sets of moments to close the equations of 
motion: the one used by Israel and Stewart |15) . and the 
one used in Ref. [TB] , and show how they lead to different 
transport coefficients. Note that, although our final equa- 
tions contain terms that were neglected in the papers by 
Israel and Stewart, see also Refs. [2"51l26j , we will still refer 
to them as IS equations. 

This paper is organized as follows. In Sec. [5] we briefly 
introduce relativistic fluid dynamics and its dynamical 
variables. The Boltzmann equation and the definitions of 
the fluid-dynamical variables from the perspective of ki- 
netic theory are introduced in Sec. [31 The orthonormal 
basis for the moment expansion and the exact equations 
for the moments are derived in Sec. 21 In Sec. [S] the 14- 
moment approximation is applied and the fluid-dynamical 
equations are derived. The choice of the moment is anal- 
ysed in Sec. |H Finally, we conclude in Sec. [7J 

Throughout this work we use natural units h = ks = 
c = 1; the metric tensor is g^ v — diag(+, — , — , — ). 

2 Relativistic fluid dynamics 

In relativistic fluid dynamics, the variables that specify the 
macroscopic state of a system are the energy-momentum 
tensor, T M ", and the particle or net-charge four-current, 
N^. Here we restrict ourselves to only one conserved par- 
ticle species or net charge. Thus, particle number and 
energy-momentum conservation imply that 

dpN" = 0, (1) 

In relativistic fluid dynamics, it is useful to define a time- 
like four- vector, u M (t,x), normalized to u^u^ — 1, and a 
projection operator orthogonal to it, 

where A» v u v = A^u^ = and A» = 3. Later on, will 
be identified as the fluid four-velocity. From now on, we 
denote the projection orthogonal to u M as A^' — A^ v A Vl 
valid for an arbitrary four- vector A^. In case of second- 
rank tensors, A 1 *", we define the orthogonal and traceless 
projection as A^ = A^ vafi A a[il where 

2 3 

Using the projection operator from Eq. ([3]), the space-time 
derivative can be decomposed as 

d„ = u^D + V„, (5) 

where the comoving time derivative is I? = u^d^, while the 
space- like gradient is V M = A^d v . For the time derivative 

we also use the notation DA — A. Applying the above 
notation, the relativistic Cauchy-Stokes decomposition is 

d ll u v =u li u v + -6A liV +cr l j tv +u llv , (6) 
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where the expansion scalar, 9, the shear tensor, , and 
the vorticity, are defined as 



(7) 

= v ( v > = i ( w + VV ) _ -eA^, (8) 



= V^u^ = i ( W - V«") . 



(9) 



The particle four-current and energy-momentum tensor 
can be decomposed with respect to as 

TV = nu^ + V, (10) 
= eu"u v - PA"" + u v W" + u"W u + tt» v , (11) 

where n = N^u^ is the particle density and e = u^T^Uy 
is the energy density. The trace of the energy-momentum 
tensor, P = — ^A^T"", denotes the isotropic pressure. 
The latter is defined as the sum of the equilibrium pressure 
Po, and the bulk viscous pressure 77, P = Po + II. The 
particle diffusion current is defined as 

V" = A»N V , (12) 

while the energy-momentum diffusion current is 

W = A^ a T a? u p . (13) 

The shear-stress tensor, , is that part of the energy- 
momentum tensor that is symmetric, traceless, and or- 
thogonal to M M , 

In local thermal equilibrium, the decompositions of Eqs. 
(jlOllip reduce to the ideal-fluid form 



K = no**". 

Tg u = e u^u u - P A^. 



(15) 
(16) 



Local thermodynamic equilibrium guarantees that the par- 
ticle density no, entropy density So, energy density So, and 
thermodynamic pressure Po, are related to the tempera- 
ture, T, and chemical potential, /_t, through an equation 
of state (EoS), i.e., Po = 7b(T, /i), from which one can 
obtain 

dP a dP 
no = ^7' S0= 9T' (17) 



and 



co = Tsq - P Q + iin 



(18) 



In general, the choice of u" is ambiguous. The frame 
where u" = u^ R — (1,0, 0, 0) is called the local rest frame 
(LRF) of matter. From the physical perspective, there are 
two natural choices which hx the LRF but at the same 
time promote it M to a dynamical quantity. According to 
the definition of Landau and Lifshitz [27] , the LRF is tied 
to the flow of energy-momentum, which leads to 



T" v u„ 



y/T^UaTtfvP ' 



and thus the energy-momentum diffusion current vanishes, 
W" = 0. 

The choice of Eckart [3B] relates the LRF to the flow 
of conserved particles as 



N" 



(20) 



which implies that the diffusion current vanishes, V" = 0. 
Sometimes it is convenient to introduce the heat flow, 



W" - h V 



(21) 



where h = (e + Po)/n is the enthalpy per particle (or per 
net charge). 

Once the four-flow of matter is specified, i.e., replac- 
ing the three independent components of W" or V" by 
u", we still have to determine 15 independent dynamical 
variables: six variables, u", e, n, and Po ; a s in the case of 
an ideal fluid, and nine variables related to dissipation, 77, 
q", and . Note that the EoS gives one additional con- 
straint and therefore reduces the number of independent 
variables to 14. 

The conservation laws (|1I2[) constitute only five equa- 
tions. Thus, to properly close the fluid-dynamical equa- 
tions it is necessary to introduce nine additional equations 
which determine the evolution of the remaining dissipa- 
tive fields, 77, q M , and ir" v . The relativistic extension of 
Navier-Stokes theory relates the dissipative quantities to 
gradients of the primary fluid-dynamical fields, 



77 = -(9, 

T 2 



(£)■ 



(22) 

(23) 
(24) 



where the bulk viscosity coefficient £, the heat-flow coeffi- 
cient Kg, and the shear viscosity coefficient r\ are positive- 
definite functions of T and \x. 

However, as mentioned in the introduction this naive 
approach leads to intrinsic problems, such as acausal sig- 
nal propagation and instabilities, and is therefore unsuit- 
able to describe relativistic fluids. The acausality problems 
were solved by introducing memory effects into the defi- 
nitions of 77, q", and 7r M ", which are no longer assumed 
to be linearly related to gradients of the primary fluid- 
dynamical variables [25PHPT, 32 , 33 , 34 , 35 . Instead, they 
become independent dynamical variables that obey dy- 
namical equations of motion (which introduce the relax- 
ation times T^, 777, and r q ) that describe their transient 
dynamics towards their respective asymptotic relativistic 
Navier-Stokes solution, 



Tn 77 + 77 = -(9 



T q q 



2j 1 a^ + ..., 



n^ 

pis 



(25) 
(26) 
(27) 



where = A£Dq a and = A^Dtt 010 and the dots 
denote possible higher-order terms. These are the type 
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of equations of motion which can also be derived from 
relativistic kinetic theory as shown by Israel and Stewart 
and others [51 fTSl[251l5B] . Causality is guaranteed, provided 
the relaxation times fulfill certain constraints [3]. 



3 The relativistic Boltzmann equation 

Let us consider a relativistic dilute gas characterized only 
by the single-particle distribution function fk = f{x^, fc M ), 
the evolution of which is given by the relativistic Boltz- 
mann equation [7], 



fc^/k = C [/] , 



(28) 



where k^ = (k , k) with k° — \/k 2 + m 2 and m being 
the mass of the particles. For the collision term C [/], we 
consider only elastic two-to-two collisions with incoming 
momenta k, k' , and outgoing momenta p,p' , 



C [f] = i J dK'dPdP'Ww 



>pp/ 



/p/p'/k/k' - /k/k'/p/p' 



(29) 



where v = 2 is a symmetry factor. The Lorentz-invariant 
phase volume is dK = gd 3 \i/ [(27r) 3 fc°] , with g being 
the number of internal degrees of freedom. The Lorentz- 
invariant transition rate Wkk/->pp/ is symmetric with re- 
spect to the exchange of the outgoing momenta, as well 
as to time reversal, 



'->pp' 



./— s-p'p 



w, 



pp/— >kk/- 



(30) 



Here, we also take into account quantum statistics and 
introduced the notation /k = 1 — a/k, where a = 1 (a = 
— 1) for fermions (bosons) and a = in the limiting case 
of classical Boltzmann-Gibbs statistics. 

The particle four-flow and the energy-momentum ten- 
sor are identified as the first and second moments of the 
single-particle distribution function, 

= (#*) , (31) 
T" v = {k^k v ) , (32) 

where we adopted the following notation for the averages 



= J dK (...) / k 



(33) 



Making use of the properties of the transition rate (j3"0"|) . 
one can show [Jj that the particle four-flow and the energy- 
momentum tensor satisfy the conservation equations (|H2p 
for any solution of the Boltzmann equation, 



9 M (k^) = J dKC [/] = 0, 



dKk v C [/] = 0. 



(34) 
(35) 



In order to identify the macroscopic variables intro- 
duced in Eqs. (fTUl ITTj) in terms of the single-particle dis- 
tribution function we decompose the momentum of the 



particles fc M into two parts: one parallel to the flow veloc- 
ity and the other orthogonal to the latter, 



k^ = E k u^ + k {f "\ 



(36) 



where we defined the energy of a particle as E k = u^k^. 
Using the above decomposition in Eqs. (|31I32[) we obtain 



T"" = (El) u^u v + -A» v {A af3 k a k p ) 



(37) 



+ u v hhJ&W ) + (E*k {v) ) + (k { >*k v) \. (38) 



Comparing these decompositions with Eqs. (jlOlllj) . we 
identify the main fluid-dynamical quantities as averages 
or moments with respect to an arbitrary solution of the 
Boltzmann equation, 

n = (E k ) , e = (El) , P + 77 = -i (A^k^k v ) , 

V = (k { ^ , W = (75 k fc <M> ) , tt^ = (k^k u) \ . (39) 

Similarly, we introduce the average with respect to the 
local equilibrium distribution function /ok, 



dK (...)/< 



Ok, 



where 



fok{x*, fc M ) = [exp (f3 E k - oq) + a] 1 



(40) 



(41) 



Although /ok satisfies detailed balance, it is not a solu- 
tion of the Boltzmann equation. The quantities ao(^ M ) 
and fi^x^) arc defined for an arbitrary non-equilibrium 
distribution function f k by the matching conditions, 

n = n Q = (£ u ) , e = e Q = (E%) . (42) 

In local equilibrium we would then identify /3q = 1/T as 
the inverse temperature and o.q = /i/T as the ratio of 
chemical potential over temperature. The matching con- 
ditions lead to 



(E k ) 5 =0, (El) s =0, 



(43) 



where (...) s = (...) — (. . .) . The matching conditions 
(|42p are convenient as they allow us to use equilibrium 
thermodynamic relations between n, e, P$, T, and /i. We 
also note that the EoS is not an additional input, but fol- 
lows from the single-particle distribution function in local 
equilibrium. 

Finally, the separation between thermodynamic pres- 
sure and bulk viscous pressure is achieved as 

Po = — (A^k^k v ) Q , 77 = - 1 (A^k^k v ) & . (44) 

We remark that (k^) Q = (E k k^) = (k^k v ~>) = 
guarantees that in local thermal equilibrium all dissipative 
currents vanish. 
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4 Moment expansion 



where the weight is defined as 



The method of moments is the most common approach 
to derive the so-called second-order theories from kinetic 
theory. In this section, we review the basic ideas of this 
method along the lines of Refs. [7J5]. Since we are in- 
terested in near-equilibrium solutions of the Boltzmann 
equation, we start by expanding /k around the local equi- 
librium distribution function /ok, 



(45) 



ui 



(21- 



— (A^k a kp) /ok/ok- (52) 



/k = /ok + Sfa = /ok 1 + /ok^k 



The coefficients olr and the normalization constants J\f^ 
can be found via Gram-Schmidt orthogonalization using 
the orthonormality condition (f5"Tj) . see Sector Ref. [T7J 
for more details. 

Finally, the single- particle distribution function can be 
expressed as 



where (f>k represents a general non-equilibrium correction. 

Following Rcf . [7] , </>k is expanded in momentum space 
with the help of the irreducible tensors 1, k^ , k^k v \ 
k^k v k x \ forming a complete and orthogonal set, 
analogous to the spherical harmonics [H] • These irreducible 
tensors are defined by using the symmetrized traceless pro- 
jections as 



oo N e 



/k = /ok i + /ok££ y&Ptr-"'^ • • • fc M 



e=o n=0 



(53) 



where we introduced the energy-dependent coefficients 



. . . jUMm) _ AHl—HmVl—Vm h, , . . L, 
i\j Fv — i—X ' yj \J\ l^m 1 



(46) 



see also Ref. [37]. The tensors k^ 1 • ••fc' im > satisfy the 
following orthogonality condition, 



dKFkk^ 1 ■■■k fMm) k^ 1 •••fc l/ ™> 
m\S r 



(2m + 1)!! 



A" 



"™ J dKF k (A a ^k a k p ) m .(47) 



Here m, n = 0, 1, 2, . . ., is an arbitrary scalar function 
of _Ek, and (2m +1)!! denotes the double factorial. 

Using these tensors as the basis of the expansion, the 
non-equilibrium correction can be written as, 



A k K <Ml ' ' ' K W>' 



£=0 



(48) 



where the index £ indicates the rank of the tensor A^ 
and £ = corresponds to the scalar A. The coefficients 

en 



«) ma y ^ e f ur ^ ner expanded in energy 7?k with an- 



l k 

other orthogonal basis of functions P kn , 



A k — L ri *kn ' 

n=Q 



(49) 



where ci^ 1 ^ are as of yet undetermined coefficients and 
Ni is the number of functions P$ considered to describe 



k 

sen to be polynomials of order n in energy, 



the t-th. rank tensor Ajf 1 ' " w . The functions are cho- 



Men — k' 



(50) 



and are constructed using the following orthonormality 
condition, 



ClKU}' /k m -r k „ - Omn, 



(51) 



H 



kri ~~Q\ / j "'ran* kri 



w _ 



and the generalized irreducible moment of (5/k, 



with 



(54) 



(55) 



(56) 



Using this notation, the expansion coefficients in Eq. 
(H^)) can be immediately determined using Eqs. flUJ and 
([51]) . For n < Ni they are given by 



^(flX—flt) 



£\ 

r=0 



(57) 



Naturally, the dissipative currents are related to the 
tensors pjf 1 ' w . According to Eqs. ([3"9")l we can identify 
them as 



pa = -377 /m 2 , 



(58) 
(59) 
(60) 
(61) 



Furthermore, the matching conditions imposed in Eq. (1431) 
can also be recast using the irreducible moments, 



P\= P2= 0. 



(62) 



The definition of the LRF corresponding to Landau's choice 
([in]) requires that 

Pt = 0, (63) 
while Eckart's definition (|20p leads to 



p%=0. 



(64) 
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4.1 General equations of motion 



while the equation for p£" is 



So far, the single-particle distribution function was ex- 
pressed in terms of the irreducible tensors p^ 1 '"^. The 
time-evolution equations for these tensors can be obtained 
directly from the Boltzmann equation by applying the 
comoving derivative to Eq. (|55|) together with the sym- 
metrized traceless projection, 



p<?i-«> = a%:::% — \ dKE^ ■ ■ -k^5f k . (65) 



Now, using the Boltzmann equation (1281) in the form 



= + 2a^a^ + [(r - 1) m 4 p 
- (2r + 3) m 2 p r + (r + 4)p r+2 ] o» v 



V-2 



rm p" r _ 1 - (r + 5)p; +1 
2 



, nv\ ■ V7<U ( 2 V) v) \ 

+ rp^ux - -V> M [m p r i 1 - p, r ^ J 

+ i[(r-l)mV- 2 -(^ + 4)pr]^ 
(2r-2)m 2 p^-(2r + 5)^^" 



(69) 



<5/k = -/ok - E^k v V v f u - E^k u V v 5h 
+ E k 1 C[f] , 



(66) 



and substituting into Eq. (|65[) . we obtain the exact equa- 
tions for p^-M . 

Since fluid dynamics does not involve tensors of rank 
higher than two, it is sufficient to derive the time-evolution 
equations for the fields p r , p£, and p£ v only. Similar equa- 
tions could also be derived for the higher-rank irreducible 
tensors, if needed. Thus, using Eqs. (|6"5j) and (|6"6"]l . the 
equation for an arbitrary scalar moment is 



Pr = C r - 1 +a<?>9+ r(f r _ x + 



D 



211 



V^_r + - ^d»W» 



+ 



1 



'20 ^20 

,2 . 



(r — 1) m p r _2 - (r + 2) p r 



3^/7 

A>0 



(r-l)pTa + 



G2r 



.//.;/ 



(67) 



Here we introduced the generalized collision term 

= a%:::% J dKE&f* ■■■wc [/] . (70) 

All derivatives of ao and /3o that appear in the above equa- 
tions were replaced using the following equations, obtained 
from the conservation laws ((I) and ([2]) , 



«o = r=— [-J 30 (M + d^) + J 20 (so + P + n)6 
U20 

+J20 (fyW - iy% - tt""^)] , 

k = tJ- [-J20 M + + J10 (£0 + Po + 71) 

D20 

+J10 (d^W - W^iifi - n^a^)} , 

u" = /V 1 (h^ao - V"A>) - — {nu» - V^TT) 

no 



(71) 
(72) 



-1 r 



"0 
n 



(73) 



where ho = (eo + -Po)/"o- The coefficients a r are functions 
of thermodynamic variables, 



Similarly, the time-evolution equation for the vector mo- 
ment is 



+ - [rmVi - (r + 3) p r+i - ia h r n] 
- iv (mV-i - Pr+i) + ^V^iT 

-A* (V X p V rX + <* h r d X -K vX ) 

+ \[{r- 1) m 2 p^_ 2 - (r + 3) ^ - 4oJfW*] 
+ 1 [(2r - 2) m 2 /£_ 2 - (2r + 3) - 5a^"] < 



+ + c^VT) < + (r - 1) p^^, 



(68) 



n 



(l-r)/ rl -/ I . --^(/ l0 G , 2 r -G 3l .), (74) 



211 



— Jr+lS — h Q 1 Jr+2,l, 

=/ r+2> i + (r-l)7 r+2 ,2, 



A) 



-./, 



r+2,1) 



where we used the notation 



(75) 
(76) 

(77) 



Im = (2 ( g+ 1) 1) ,, J dKE^T 2q (A^k a kp) q fo^ (78) 



" 9 (2? + l)!! 

C nm = JnoJmO ~ Jn— l,oJm+X»0 
&nq — Jn-\-l,qJn — l,q i^Jnq) 



M^ fc «fc/3) /ok/ok, (79) 



(80) 
(81) 
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Thus, we have obtained an infinite set of coupled equa- 
tions containing all moments of the distribution function. 
Note that the derivation of these equations is independent 
of the form of the expansion we introduced in the previous 
subsection. 



5 The 14-moment approximation 

In order to obtain the macroscopic equations of motion in 
terms of the fluid-dynamical variables that appear in the 
particle four-current and energy-momentum tensor, the 
generic hierarchy of the coupled moment equations must 
be truncated. To this end, Israel and Stewart made the 
so-called 14-moment approximation: they truncated the 
expansion of the distribution function and matched the 
non-equilibrium corrections to the dissipative currents 77, 
V 1 , W», and tt^. 

In this section, we show how the 14-moment approxi- 
mation emerges from the general moment expansion pre- 
sented in Sec. SJ First, we neglect irreducible tensor mo- 
ments of rank higher than two, i.e., p^ 1 '"^ — for £ > 3 
in Eqs. (|68I69|) . Such irreducible moments cannot be con- 
structed purely from first-order gradients of equilibrium 
fields [T7]. This means that they lead to terms that are 
of higher order in gradients or contain higher powers of 
dissipative quantities in the equations of motion. 

Next, in the expansion (I53|) of the distribution func- 
tion we include the first three scalar moments, namely 
p = —3II/m 2 , p\ = 0, and p 2 = (the last two scalar 
moments vanish due to the matching condition, but must 
be included since they were used to define ckq and /3q), the 
first two vector moments, Pq = and Pi = W* 1 , and the 
first second-rank tensor moment p^ — tt^ . This implies 
that No = 2, Ni = 1, and N 2 — 0, while all other mo- 
ments appearing in the expansion are dropped. Choosing 
either the Eckart or the Landau frame, we can eliminate 
either V 1 or W^, respectively. Let us note that, so far, this 
approach is completely equivalent to the matching proce- 
dure of Israel and Stewart. The ambiguity in the transport 
coefficients emerges only at a later stage, when choosing 
moments of the Boltzmann equation to supply the equa- 
tions of motion for the dissipative currents. 

The restriction to the aforementioned moments affects 



that Ak, A 



£>, and A<T> 



from Eq. (fi9")) can be expressed 



solely in terms of II, V 11 , , and 7r 



No 



A k - E CnPS - col® + clPM + c 2 P$, (82) 

n=0 
Ni 

# ) = E# )p £!* 1 # >P & ) + # >F ®> (83) 

n=0 
N 2 

\ (/*«) _ „{p.u) p(2) ^ Uiv) p (2) 



(84) 



n=0 



where the tensors 



are given by Eq. (|57|) . while 



those which do not appear in the above equations are set 



to zero. According to Eq. (|57[) the scalars cq, c\, and c 2 
are proportional to the bulk viscous pressure, 



CO = 


HI 


(85) 


Cl = 




(86) 


C2 = 




(87) 



The vectors Cq and cf^ 1 are given by a linear combination 
of particle and energy-momentum diffusion currents, 

^^"W (88) 
c<"> = V^aftW™ + W»a^N {1) , (89) 



while is proportional to the shear-stress tensor, 



■™ J 2 ) 



Let us recall Eq. ([5U| and for any £ > we set 



McO — u 00 " L ' 



while 



(90) 

(91) 

(92) 
(93) 
(94) 



The orthonormality condition (|5ip implies that the nor- 
malization constant is 



p(0) _ J0) p v 
P (i) _ M) F , , „(i) 

p(0) _ J0) F 2 , J0) p ,10) 
Mc2 — a 22 ^k + a 21 ^k + a 20 . 



,(0) 



and 



,(0) 



'ii 

,(0) 



,(o) D 



i. 22 



■^10 
^00 


, (oi?) 


2 


G12 


fl (0) 

"20 


£>20 




(0) 
"22 


Dio 




•^00 


Dio 



72 



no 
'ii 



J20D20 + J30G12 + J40-D10 
■hi 



^ (off) 

21 ^ ' 



2 72 

_ J 21 



Ail 



(95) 

(96) 

(97) 
(98) 
(99) 



Furthermore, using the orthogonality relation (|47p together 
with Eqs. (|48H50[) one can easily show that 



Ne n 

f--^ = ^EE c « w ' 

n— m— 



a nm Jr+m+U,l- 



(100) 



Applying the truncation scheme required by the 14-moment 
approximation we obtain that all scalar moments, p r , be- 
come proportional to the bulk viscous pressure 77, 



No n 

Pr = y ] y ] ^v" 1,111 
n=0 m=Q 



fl (0) J j. n 



(101) 
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Similarly, all vector moments, p%, are proportional to a 
linear combination of and W* , 

Nx n 

rf^EE c<f>a« J r+m+2)1 = yYv'+iyW, (102) 

n— m— 

and, finally, p£ v is proportional to tt^ , 

N 2 n 



'^EE c ™" y>a «™ = 7> M ^ (103) 



n— ?n— 



Now, using the previously obtained results we prove 



that 



where 





if 1 


= A n J r o + B 




CnJr+2,0, 


(104) 






= AyJr+2.1 d 


- By Jr+3.1 




(105) 




= AwJr+2,1 ~ 


h Bw^r+3,1) 


(106) 






= 2„4 7r J r+ 4 j 2, 






(107) 


An 




3 


^30 




(108) 




m 2 J 2 qD 20 4 


■ ^30^12 + 




B n 




3 


G23 




(109) 




m 2 J20-D20 ■+ 


■ J30G12 + 


J40D10 ' 


C n 




3 


D 20 




(110) 




m 2 J20D20 d 


- J30G12 + 


JaqDxq ' 


A v 




J41 

— — , «4w — - 
^31 


J31 

D 31 ' 




(111) 


By 




^31 R 
-75—, % - 


J21 
Dai 




(112) 






1 

2 J42 






(113) 



We remark that, since the matching conditions were al- 
ready imposed, one can prove that 7^ = 7^ = 0. 

Equations (|101|) . (|102|) . and (|103j) are the main result of 
the 14-moment approximation. Such relations guarantee 
that any irreducible moment of the distribution function 
can be expressed in terms of the dissipative currents ap- 
pearing in and T^ v . This is also what Israel and Stew- 
art achieved by their matching procedure. Consequently, a 
closed set of fluid-dynamical equations can always be de- 
rived. This happens because the reduction of dynamical 
variables was done directly in the single-particle distribu- 
tion. On the other hand, this truncation also leads to an 
ambiguity in the derivation of fluid-dynamical equations 
since, for example, the equation of motion for the bulk 
viscous pressure can be obtained from p r for any r. We 
will come back to this point in Sec. [5] 



5.1 The collision term 

In order to express the collision term (1701) in terms of 
the fundamental fluid variables, C[f] is linearized in de- 
viations from the equilibrium distribution function. Sub- 
stituting Eq. (j4"S"]) into the linearized collision term, we 



obtain 



C£i" w - - / dKdK'dPdP'W kkl ^ ppl f Qk f okl f 0p f 0pl 



x E^k^ ■ ■ • fc« (0 p 4- V - <f>* - <t>v) , (H4) 

where the cf>'s are given in Eq. (|4"5)l . In order to specify 

CV_i, Cy™\, and C^°fi for the general equations of motion 
(|67l l68l l69|) we project the collision term as 

CV-i = it Ml • • • u^C r j^_* , (115) 

cl a \=Al l u^---u llt C^> H , (116) 

In the 14-moment approximation we start by substi- 



tuting Eqs. dH]) and (@9J) into Eq. (1TT4[) and obtain 

CV_i = Cn x nX r -a,i, (118) 
C^ a = By^X r _ 2 , 3 + BwWXr-ifr (119) 



a 



r— 1 



(120) 



Here X r i, A r .3 and X r! 4 are coefficients of the following 
rank-4 collision tensor, 



X 



1 



- / rf^rfif'rfPdP'W kk ^ pp/ / ok /ok//op/op/ 



x EtfPW (p a p p + p' a p' -k a k l) - k' a k' p ) 



(121) 

This collision tensor is symmetric upon the interchange of 
indices (/i, v) and (a, /3), and also traceless for the latter 
indices, 

x „vaf) = x (jH>H<xfl) t X? vaf> g af) = 0. (122) 

These properties lead to a spatially isotropic tensor con- 
structed using the four- velocity u 1 *, the transverse projec- 
tion A^ , and different scalar coefficients X r ^ as 

X p va fi = ( XrA u»u v + X ri2 A" v ) (u a u p - ^A afj 

+ AX r ^ A v ^ {a +X rA A^ va ^ (123) 

where the coefficients of the above decomposition are ob- 
tained from the following contractions, 



X r<1 = Xr a0 u^u a up 

X r ,2 = ^Xjt Vaf, A llv U a Up 

X r , 3 = -X^u^A^up 
X r ,4 — -^Xl^ val3 A^ va p . 



(124) 
(125) 

(126) 

(127) 



The evaluation of the coefficients X r ^ requires the de- 
tailed knowledge of the differential cross section. As an 
example, these functions are evaluated for a classical gas 
of massless particles with constant cross section in Ap- 
pendix 121 
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5.2 Equations of motion 

In order to close the conservation laws Eqs. §T§ and (J2J 
we need additional equations for the dissipative currents 
which can be obtained from the exact equations of motion, 

Eqs. dsn Ell ESI)- 

For any r > 0, Eq. (|67|l for the scalar moment leads 
to an equation of motion for the bulk viscous pressure. 
Replacing p r = 7^77 according to Eq. (|101[) in Eq. (joTl) 
and collecting all terms we obtain the so-called relaxation 
equation of the bulk viscous pressure, 



c 



•n 



5 r nn m 



(128) 

Here, we introduced the relaxation time of the bulk viscous 
pressure, rjj, and the bulk viscosity coefficient £ r as 



-3,1 



r ai 0) 



(129) 



where Cn was defined in Eq. (I110[) . Similarly, 7^ was de- 
fined in Eq. (fT04|) . X r<1 in Eq. ([I2g]> . while a[ 0) was given 
in Eq. (|73|) . The other transport coefficients in Eq. (|128l) 
are defined as 



<nw 



' nv 



1 

w 

1 



(r- 



G*2r 



Po—^TB ~R — 



(r -1)7,1 



■A) 



£>2 



0A) 



A w — 



1? 





,_lC7r-l 


da Q 


f dp 






dao 


f a/? 



pr 

l nw 



pr 



J nn 



11 



7, 



1? 



1? 

HqD 



w 



G-2r 



(r- 



1)7?_ 2 + ^ 
-^20 



20 



11 



>30 



57. 



'20 



3 g7j 
da 

di n ~ 



'10 



97 



J3 



0A> 



0A> 

(r + 2) 7^-3^1 



37" 



(130) 
(131) 
(132) 
(133) 
(134) 
(135) 
(136) 

h 

(137) 



Note that these coefficients are independent of the colli- 
sion integral. We also point out that here we follow the 
notation of Rcfs. [25,22] for the coefficients. Furthermore 
the choice of the LRF eliminates terms involving either 
(for Eckart's choice) or W 1 * (for Landau's choice). 



In the very same manner we get a relaxation equa- 
tion for both the particle diffusion current and the energy- 
momentum diffusion current. This equation follows from 
Eq. ([SSI) using pjf = t^V" + -y^W, where rank-3 ten- 
sors p^ vX for any r > were neglected. Thus, after some 
calculations we obtain 



v 

1 1, 



v. 



w 



+ 



•v 



' B 2 h 2 



Va 



- ^XwwWvO^ + \ r vv V v a^ 
-i>™5 r ww W»6 + 6 r vv V»6 

- T r qn IIu» - T r q ^ v u v 

+ A^TTV^o + A^Tr^V.ao, 



(138) 



where we defined the relaxation time of the particle diffu- 
sion current, t v , of the energy-momentum diffusion cur- 
rent, Tm, and the heat conductivity coefficient k£, 



Mr 



7, y 



W 

It 



a' 



a {1) 
r v~~v~ n oPo- 



Furthermore, 



W 



7r v 



w 

Ir 



al 



(139) 
(140) 
(141) 

(142) 



and By, Bw were defined in Eqs. (|112l) . such that By = 
-h B w , while 7^, 7^ were defined in Eqs. (jl05ll06j) . In 

addition, X r ^ was defined in Eq. (|126p . while a^ and 
a£ were given in Eqs. ([75)1 and ([77]). The two relaxation 
times are not independent but are related to each other as 
Tyy = —hoip^Ty. The relaxation equation is written such 
that it is straightforward to rewrite it in either the Eckart 
or the Landau frame. Moreover, using the definition (f2"l"|) 
of the heat flow it is clear that the above equation is a 
relaxation equation for this quantity. 

The coefficients in Eq. (|138[) which only exist in either 
the Eckart or Landau frame are 



A 



WW 



'WW 



* [m a (2r- 2)7^2 



57 r K i>: 

(2r + 3) 7 f 



J: 



dl7 



J: 



20 



- 5<] 
dao 



w 



dao 9/3o 
(r + 3) 7 r-4^], 



1 



(143) 

[m 2 (r-l) 7r W 2 
(144) 
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and 



1 



X vv = TTv h 2 ( 2r - 2 ) ~fr-2 - (2r + 3) 7r y ] , 

°7r 



(145) 



°vv — 



71qD. 



-i r 



'30 " 



X 



3 7 7 



J20 
J20 



v 



da 



dj3 



1 



+ h 2 (r ~ 1) 1Y-2V 1 - (r + 3) 7l 



r ] ) 



(146) 



while terms and coefficients which are not affected by ei- 
ther choice of the LRF are 





1 


T r = 
' q-n 


"tT 




1 


T q n — 


~3t7 



(r-i) 7 ;_ 1 + A 



1 a/3 



o 2n dlr-1 o 9 7^+i n 

-3a r + m p — — Po 



dpo 

1 



c g 77 



7 r V 
1 

"7 



7, 



(mV-t-7^) 



7r ^ \, da + 0/3 o 



1 

3t7 



/^-i , ,-1^-1 



/ ^7^1 h -l d lr+l 

\ da + "° 9/3 



(147) 

(148) 
(149) 
(150) 
(151) 

(152) 



The relaxation equation of the shear-stress tensor fol- 
lows from Eq. (|69|) by replacing pf v ^ = 7 ^7r M ", and ne- 
glecting rank-3 and rank-4 tensors, p^ vX — and p^ LuXK = 
0, 



2rf 



- 2\l^u v ^ a + 2^uo^ a - 2S r ^ v 6, (153) 



where we defined the relaxation time t£ for the shear- 
stress tensor and the shear viscosity coefficient rf, 



7; 



(2) 

rf=<^. (154) 

7r 



coefficients in Eq. (|153[) are 
1 



Kn — 



-[{r-l)m i 1 ^ 2 -{2r + i)m 2 1 n 



15 7 , 

+(r + 4)7& a ], 
1 



5 7 ? 

(r + 5) 7 £, 
1 



9-y w 

V7 a u 7r+l 



'ttV 



d/3 



2 . \ 2 j ^7r- 1 



57? 

(r + 5) 7 ^ +1 -/3 
1 



pr 



a ttW 



57 
1 

~H 
1 

5^ 



2 Vf w \ 
- (m 7r _! - 7r+1 ) , 



12V V \ 

- (m 7r _! - 7r+1 ) , 



9«o ° 0A> 



Kv — 



1 

57? 



^7r-l , , -l djr-1 

da d/3 

V 



"0 



<9a 



(155) 
(156) 

(157) 
(158) 
(159) 

(160) 
(161) 



K« = [(2r - 2) m 2 7r % - (2r + 5) 7?] , (162) 



2 7 - 



^7r" 



57? 

a 7 ; 



7 LI T 

■i[(r-l)m^ 2 -(r + 4) 7r '] 



(163) 



In the above relaxation equations we have expressed 
the proper-time derivative and spatial derivative of the co- 
efficients from Eqs. ([H? H WI \ ITO ^ using Eqs. ([7T] [7l2l [75]) . 
Therefore, for any coefficient 7? , 7) T, 7^, or 7 ^, collec- 
tively denoted by ip (ao, j3o), we used the following formula 
for the proper time derivative 



^~ £20 



6V 3^' 
J20T; 1- ^iottj- «o 



J3. 



9^ 



J 2 



9-0 



9a a/? , 

while for the spatial gradient of ip we used 



(164) 



(165) 



Here, A* was given in Eq. (|113p . 7 ^ in Eq. (|107l) . X r ^ 4 in 
Eq. (T1271) . while a[ 2) was quoted in Eq. ([75|. The other 



Note that these two equations follow from the equations of 
ideal fluid dynamics and we neglected terms from the gen- 
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eral conservation equations proportional to the dissipative 
fields or their derivatives. 



6.2 Equations of motion directly from the dissipative 
currents 



6 Choice of moment and coefficients in the 
massless limit 

As was shown in the previous section, once the 14-moment 
approximation is applied, any moment of the Boltzmann 
equation will lead to a closed set of equations which looks 
formally the same, see Eqs. (fT28)l . (fH8]l . and (TT53J). This 
is immediately apparent in these equations by the explicit 
dependence of the transport coefficients on the index r. 
Barring any miraculous cancellation, already at this point 
it is obvious that their values will in general be different 
for different r. Thus, the 14-moment approximation leads 
to an ambiguity, since it gives rise to an infinite set of 
equations to describe a finite set of macroscopic variables. 

In order to make this clear, we shall consider two dif- 
ferent choices for the moments of the Boltzmann equation: 
the first one is the traditional choice of Israel and Stewart 
[15] and the other one was recently proposed by Denicol, 
Koide, and Rischke (DKR) jTB] , following Grad's original 
idea. 



6.1 Equations of motion of Israel and Stewart 

Israel and Stewart assumed that the equations of motion 
for the dissipative currents could be extracted from the 
second moment of the Boltzmann equation [SHOUTS]. 

d^(k^k v k x ) = J dKk"k x C[f], (166) 

with the equations for 27, V M or W 1 *, and tt 1 *" obtained 
using the following projections of the above equation, 

u v uxd„ (k"k v k x ) = u„u x [ dKk v k x C [f] , (167) 



A a x u u d^ (k»k»k x ) 



Af x d,(k^k x ) = Af x 



dKk v k x C[f], (168) 
dKk u k x C[f], (169) 



respectively, together with the 14-moment approximation. 

As a matter of fact, Eqs. fl67|) . ([1081 and ([TBS]) can be 
identified as the equations for p3, and p\ . These equa- 
tions have already been calculated with the 14-moment 
approximation, Eqs. (|128[) . (|138f) . and (|153l) . Thus, using 
the indices r = 3 (for the scalar moment), r = 2 (for the 
vector), and r = 1 (for the irreducible second-rank tensor), 
we obtain the IS equations. 

Even though this choice of moments was never clearly 
justified, this prescription is widely employed in relativis- 
tic kinetic theory. However, it was recently found that, 
at least for some cases, the IS equations are not in good 
agreement with the numerical solution of the Boltzmann 
equation ri8lfT9l[2^[2T1 [22 l [23 l . Also, the transport coeffi- 
cients obtained by Israel and Stewart do not coincide with 
quantum- field theoretical calculations [23] . 



In the second choice, the equations of motion for the dis- 
sipative currents are obtained from the moments po, pg, 
and Pq U which exactly correspond to the dissipative cur- 
rents, see Eqs. (|58]l. (|59]l. and ([6"T]) . Here, we have already 
fixed the LRF in accordance with the Landau picture. 
Then, the equations of motion for the dissipative currents 
emerge directly from their definitions as 



n=-\m 2 J dKSf k , 
= / dKk w Sf k , 



*■</*"> = / dKk { »k v) 5h 



(170) 
(171) 
(172) 



As already mentioned, these equations are related to the 
equations for p = — 377/rn 2 , p£ = V\ and p^ v = ir^" , 
that is, Eqs. ([67J). ([55]). and flOO]), for r = 0. In this sce- 
nario, the fluid-dynamical equations (11I2[) are closed by 
Eqs. ([T28"|) . fllUg]) . and ([T55|) , with r = 0, which correspond 
to Eqs. (11701 . (1TTT]) . and ([172]) . once the 14-moment ap- 
proximation is applied. 

It was shown that this method can successfully repro- 
duce the numerical solution of the Boltzmann equation 
for the simple one-dimensional scaling expansion [16) . It 
is also important to mention that the transport coeffi- 
cients of this kinetic calculation are consistent with those 
calculated from quantum field theory with the method 
proposed in Ref. [21] , 



6.3 Comparison of choices 

In order to understand the difference between the two 
approaches discussed above, we calculate the coefficients 
Ph = -C/rh, & = V r K, and ft = K r q / foffih®. 
These coefficients, normalized by the pressure or parti- 
cle density, are shown in Figs. [JJ [2] and [21 respectively. 
The calculations were done for a classical gas with fixed 
chemical potential, p = 0. 

We see that both calculations converge at low tempera- 
tures but deviate considerably at high temperatures. This 
behavior should be qualitatively the same for any choice 
of moment because all irreducible moments of the same 
rank converge to the same values in the non-relativistic 
limit (multiplied by a different power of the mass). Thus, 
differences between the choice of moment will only appear 
in the relativistic limit. 

The coefficients in the ultrarelativistic limit, m/T = 
0, for a classical gas with constant cross section, can be 
calculated analytically. These are collected for the shear 
viscosity and particle diffusion in Tables [TJ [2] and [3] Note 
that, in this limit, the bulk viscous pressure vanishes and 
was not considered. For the relaxation times, Ty and t w , 
and transport coefficients, r\ and Ky, the differences are 
of the order of 10 — 20 %, but for other coefficients the 
differences can be more significant. 
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Table 1. The coefficients for the shear-stress tensor in the two approaches for the classical gas with constant cross section in 
the ultrarelativistic limit. 
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Table 2. The coefficients for the particle and energy diffusion in the two approaches for a classical gas with constant cross 
section in the ultrarelativistic limit. 
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Table 3. The coefficients which couple shear stress and particle or energy diffusion in the two approaches for a classical gas 
with constant cross section in the ultrarelativistic limit. 
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Fig. 1. The coefficient At normalized by the pressure Po- The 
cases r = 3 (dashed line) and r — (solid) line correspond to 
the choices by Israel and Stewart and by Denicol, Koide and 
Rischke. 
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2. The coefficient Ar normalized by the pressure Po- The 
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the choices by Israel and Stewart and by Denicol, Koide and 
Rischke. 



The coefficients AL.,.^., and r q7r in Eq. (|138D are cal- 
culated both in the Eckart and Landau frames, see Table 
3. For example, in the Landau frame we only have an 
equation for and hence Ay w = Xg„ is given exactly 
by Eq. ()151[) . while in the Eckart frame we only obtain 
an equation for W so that, X 7 kW — A^/V'w- However, if 
we use the definition of the heat flow in either frame, i.e., 

= —hoV* 1 or q M = W M , then these coefficients lead to 
the same values for a classical gas where ip^ = —h^ 1 . 



7 Conclusions 

In this work we have reviewed the 14-moment approxi- 
mation proposed by Israel and Stewart and discussed the 
ambiguities of this approach. We started by introducing a 



general expansion of the single-particle distribution func- 
tion in terms of its moments. For this purpose, we con- 
structed an orthonormal expansion basis which allowed us 
to establish exact relations between the expansion param- 
eters and the moments of the distribution function. We 
then proceeded to derive the exact equations of motion 
for these moments. 

Next, we showed how the 14-moment approximation 
can be obtained as a truncation of this general expansion 
of the distribution function. We proved that, once the 14- 
moment approximation has been applied, it is possible to 
derive an infinite number of fluid-dynamical equations, all 
having the same general structure but with different trans- 
port coefficients. This means that the 14-moment approx- 
imation is not able to provide a unique theory of fluid 
dynamics and, in this sense, is ambiguous. In Sec. [B] we 
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Fig. 3. The coefficient (3y normalized by the pressure Po- The 
cases r = 2 (dashed line) and r — (solid) line correspond to 
the choices by Israel and Stewart and by Denicol, Koide and 
Rischke. 



analysed two different choices for the moment equations: 
the one corresponding to Israel and Stewart |15j . and the 
other one to that of Denicol, Koide, and Rischke |16j . It is 
also worth to mention that in this derivation we obtained 
terms that were neglected in the original work of IS [15 , 
as was already presented in Ref. [251125] . 

We also remark that the solutions of the IS equa- 
tions were already compared to the numerical solutions of 
the Boltzmann equation for the so-called Bjorken-scaling 
problem in Refs. |181[T9"1[20"] and for the relativistic Rie- 
mann problem in Ref. [2TU22H25j . It was demonstrated 
that IS theory is in relatively good agreement with the 
numerical solutions of the Boltzmann equation only if 
the Knudsen number is sufficiently small. Note that these 
comparisons did not include all non-linear terms and trans- 
port coefficients derived in this work. On the other hand, 
in Ref. [16] it was shown that, in contrast to IS theory, 
the direct method gives a much better agreement with 
the numerical solution of the Boltzmann equation up to 
very large Knudsen numbers. 

Before closing we mention that recently the method 
presented in this work was extended to include 14 + 9 x n 
moments. It was explicitly shown how to successively im- 
prove the expression for the transport coefficients by ex- 
tending the number of moments from n = to n = 1,2, 
and 3 |17j . Furthermore, it was also shown that the equa- 
tions of motion can be closed in terms of 14 dynamical 
variables without making use of the direct truncation of 
the moment expansion, the 14-moment approximation. 
This was obtained by a separation of the microscopic time 
scales and a power-counting scheme in Knudsen and in- 
verse Reynolds number. The equations of motion can be 
closed in terms of only 14 dynamical variables, as long as 
we only keep terms of second order in Knudsen and/or 
inverse Reynolds number. 



This work was supported by the Helmholtz International Cen- 
ter for FAIR within the framework of the LOEWE program 
launched by the State of Hesse. The work of H.N. was sup- 



ported by the Extreme Matter Institute (EMMI) and the Aca- 
demy of Finland, Project No. 133005, that of E.M. by the Hun- 
garian National Development Agency OTKA/NFU 81655. 



A The collision integral in the massless limit 

In this Appendix, we calculate the collision tensor defined 
in Eq. (|121j) . For a classical gas with constant cross section, 



X 



i / dKdK'dPdP'W^ ppl f^f 0W 



x EWtf (p a p p + p' a p' p -k a k l) - k' a k lfi ) 



(173) 



First we define the total cross section as 



a T (s) = - 
v 



2tt sin <9 s d<9 s a (s,6> s ), (174) 



where a (s, S ) is the differential cross section, s is a colli- 
sion invariant, i.e., a Mandelstam variable, and O s is the 
scattering angle, 



s = (k» + k'rf = (p» + p'»f , 



& s = arccos 



(^-fc^)(p M -f/ M ) 
(k^ - fc'^) 2 



(175) 
(176) 



The transition rate Wkk/-»pp/ is written in terms of the 
differential cross section as 

Wkk/^pp/ = (2tt) 6 s a (s, 9.) 6* (k» + k'» - p» - p'») . 

(177) 

In order to simplify the calculations we divide X^ ua ^ 
into gain and loss parts, 



(178) 



where 



Gr aP = ~ / dKdK'dPdP'W^ ppl f ok f v 

x E^k^k" (p a p p + p' a p' p ) , (179) 

c ^ap = I f dK dK'dPdP'W^ ppl foKfov 

x Elk"k u {k a k p + k' a k' ) . (180) 

The tensor C^ va ^ can be directly integrated and written 
in terms of the total cross section, 

qmP = 1 j dKdK ' fmfok , UT ( s ) ^ s {s _ 4m 2 ) 

x Elk"k u (k a k p + k' a k' ) , (181) 
where we used 

1 f d 3 p d 3 p' 



v J p" p' 
1 



„ „) sa(s,O s )Si k ,^ pp> 



or (s) y/ s (s — Am 2 ). 



(182) 
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For the tensor Q^ va ^ we first introduce the total momen- 
tum and corresponding projection orthogonal to it, 



P% = k^ + k'v = p>* + p'f, 



1 rji r rj~ 

S 



(183) 
(184) 



Now, the p— dependent part of the integral can be written 

as, 



1 f d 3 p d 3 p' 



v J p° p'° 



■ s a 



(s,O s )p a P ^t k ^ pp , 



B 1 P^P^+B 2 A^, 



(185) 



where in the center-of-momentum frame in which P!^ 
(•v/i, 0,0,0) and A^,p M p v = - |p| 2 , we obtain 



B i = (s) A/s(s-4m 2 ), 
B 2 = -^a T (s) (y/~s - 4m 2 



(186) 
(187) 



In the massless limit, the above results simplify con- 
siderably, 

L ^ a! 3 = 1 J dKdK ' fmfm , S(TT ( s ) 

x EWW (k a kP + k' a k ,f) ) , (188) 



^ J dKdK'fokfok' s a T (s) 
WW (P$P? - S -g afi ) . (189) 



From here on, we will consider only the case of constant 
cross section. Then, using s = 2 (fc M fc^) = 2 (p^p 1 ^), we 
directly obtain 

0? af> = a T (E r k k^k^k a k^k K ) (k K ) 

+ ( 7 T ( J B^fc^ K ) <fc Q fc^ K ) , (190) 



and 



2 

— < 

3 

4 

+ -a T 
2 

+ ttT 



[E-^k^k k )q (k fc^fc K ) 



- -erg"? (E^k v k K k x ) Q (k K k x ) Q , (191) 

Finally, using the definition of the thermodynamic inte- 
grals from Eq. (|78|l we obtain 

A r — g 1 r+5 Jl > K ' g i r+4 1 2,k 

U T j^vnjap a T a jIJ,vk\ r 
~ "g* i r+3 i 3,« - -g - ^ V+4 J 2,kA- (,192J 



Therefore, the different projections are given as 

X r .i = X^ af} U^U u U a Ui3, 



-<JT 



{Ir+5,()IlO + ^r+3,0^3o) _ 8/ r +4,i/21 



(193) 



X r ,3 = ^X^ ual3 u^A ua u f3 



CT'T 

— [J r+ 5,l/l0 - 4/ r+4 ,l/ 2 l - ir+3,l-f3l] , (194) 



and 



2er T 



[-^■+5,2^10 + 4/ r +4,2^2l] 



(195) 



In order to calculate the coefficients in the massless 
Boltzmann limit, we use the following formula for the ther- 
modynamic integrals 



L n+r,q 



(a ,j3 ,m -> 0) 



P (r + n + l)! 
2/3 r+ "- 2 (2q + l)!!' 



where Pq = ge a °(3 4 /7r 2 , hence 



A r ,3 — 



a T P'i(r + A)\ 



(r 2 + 3r + 2) 



a r P 2 (r + 4)! 
18/?S+ 2 



(r 2 + 7r + 6) 



^=- ^£ + t 5)! (^ + io). 



45^ 



(196) 
(197) 
(198) 
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